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Abstract 

We introduce basic aspects of new operator method, which is very suit- 
able for practical solving differential equations of various types. The main 
advantage of the method is revealed in opportunity to find compact exact 
operator solutions of the equations and then to transform them to more 
convenient form with help of developed family of operator identities. On 
example of non-linear first-order DEs we analyse analytical and algorith- 
mical possibilities for solutions obtaining. Different forms of solutions for 
first-order DEs are given, including for some integro-differential equations 
and equations with variational derivatives. We describe new algorithms 
for direct computing the solutions with help of computer algebra system 
(CAS). We also discuss recipe for finding new solvability conditions, which 
allow to enlarge DE solving abilities of existent CAS. 



1 Introduction 

In recent years it becomes obvious a reanimation of intense interest to further 
development of exact methods for solving differential equations. The main rea- 
sons here are not only in expansion of science areas where differential equations 
play prominent role, but also in impetuous development of computer algebra 
systems (see [l]-[3]). 

There are several universal approaches to solving DE. The most promising of 
them are Lie symmetry methods (see [4] ) and advanced methods for the linear 
equations which are based on differential Galois theory [5] with the first usable 
algorithm developed by Kovacic [6] , [7] . But in computer algebra very often un- 
derlying mathematics is far from the methods normally taught to mathematics 
students for the same problem [1]. 

The aim of this contribution is to show on examples of first-order equations 
the possibilities for solving differential equations with new method based on 
Dyson's operator solution representation in form with time-ordered exponent. 
This method allows to find instantly the exact operator solution if we succeed 
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in reformulating given problem to the first-order linear one. This solution is 
formal enough. The next step of the approach is directed to transformation 
the operator solution to combinations of some operators with relatively simple 
properties, which would permit to carry out all the operators and obtain as a 
result the solution of equation " by quadratures" or to more practical, calculable 
or usual (standard) expression. 

The main advantages of the method are in its compactness, clarity and 
simpleness. Since we often can obtain some form of exact solution and due to 
the fact that it is much easier to study the "prepared frog". It is essential as 
well, that the operator method is algorithmic. 

It is clear that there are some correlations between operator method and oth- 
ers common methods. But we believe that in introductory paper it is untimely 
to analyse such correlations. 

We choose the first-order equation as a testing area for the operator method 
because it is most important differential equation in many branches of science 
and because main aspects of the method are quite easily generalized for other 
cases. 

There are some known solvability conditions for the first-order equations (see 
for example [8]). But in the situation, when we are far from be able to solve 
any equation, the pursuit of new such conditions is very important. 

To avoid induced complications we choose a nonrigourous manner of expo- 
sition inasmuch we always can verify obtained conclusions by substitution to 
initial equation. Analytical restrictions are obvious from context as a rule. 

The paper is organized as follows. In Sec. 2 we introduce basic ideas of 
the operator method. The applicability of the operator method for non-linear 
problems is discussed in Sec. 3 on example of first-order non-linear ODE. The aim 
of this section is to introduce on particular case what we will do for the general 
cases. In Sec. 4 we consider an idea of CAS solving algorithm, which is based 
on operator solution. We believe that this type of direct solving algorithms 
can play a noticeable part in future CAS developments. The main tools of 
operator method: the transformation identities are briefly described in Sec. 5 
and important properties of generalized shift operators are considered in Sec. 6. 
In the Sec. 7 we demonstrate algebraical way to simplify operator solutions. Such 
type of solutions can found more efficient direct CAS algorithms. The expression 
of the general solution of the equation through particular solution is considered 
in Sec. 8. This section carries a methodical load too. Here we demonstrate 
the possibilities of operator method in obtaining solutions for more complicated 
equations. It is shown in Sec. 9 that the solution of DE can be expressed in terms 
of arbitrary function. Among another consequences it brings some recipes for 
finding of solvability conditions that allow to enlarge DE solving abilities of 
existent CAS. 



2 



2 Base of the operator method 



The main point of the operator method is based on solution of the first-order 
linear equation obtained by Dyson in the form with time-ordered exponent. 

For the first- order equation with respect to t (in ordinary or partial deriva- 
tives) for u(t, p) : 

^Pl=L(t,p)u(t,p) + f(t,p) (1) 
with the initial condition 

u(a,p)=v{p), (2) 

where L is a linear operator (independent variable p is one dimensional or many 
dimensional), the solution is given by (t > a) 

u(t, p)=f cxp{ f dr L(t, p)}v(p) + f \ f exp{ f d( L((, p)}f(£, p) . (3) 

J a J a J £ 

Here the operator T is a "chronological" (ordering) Dyson's operator, repre- 
senting a product, ordered (in order of unincrease) along the co-ordinate t , of 
the operators, which appear in the integrands in the expansion (see below) of 
the exponent in an operator series: 

f[L(T 1 ,p)...L(T m ,p)} =L(T ai ,p)L(T a2 ,p)...L(T am ,p) . 

T a\ ^ T ct2 — ' ' ' — ''dm 

This chronological operator has to be introduced because the operators in the 
integrands with different values of the variable r may not commute. The solution 
of the problem (1), (2) for t < a is similar to (3), but instead of the operator T 
it is necessary to place the operator T , which represents a product, ordered 
in reverse order. 

We can define chronological operator exponent by another way, e.g., through 
its differential properties or by way of expansion in an operator series. We can 
find that 

f exp{ f dr L(r, p)} = 1 + [ dr L(r, p) f exp{ f d£ L(£, p)} , 
were f is the identity operator. Father iterations lead to 
f exp{ f dTL{T,p)} = l+ I dTL(r,p) + ... 

J a J a 

h / dn / dr 2 --- dr n L(n, p) L(t 2 , p) . . . L(r n , p) + . . . (4) 

J a J a J a 
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The most important here is the fact that (3) is really the exact solution of 
the equation (1) with the initial condition (2). One can very easily verify it by 
means of differentiation of (3) with respect to t . 

It is obvious that 

t exp{ f dr L{t, p)} f exp{ f dr L(t, p)} = 1 
and, e.g., when a < t < b 

f-b f-t f-t 

fexp{/ dr L(t, p}} f exp{ / dr L(t, p)} = f cxp{ / dTL{r,p)}. 

J a J b J a 

It must be stressed here that we do not involve any premises about nature of 
u(t, p) , it can be both scalar and vectorial function, and it can depend on any 
number of independent arguments. For operator L (and its "filler") as well - 
the only fundamental requirement - it must be a linear operator. Very often 
non-linear problem is reduced to linear differential equation. If we succeed in 
reformulating given problem to first-order linear one, we then instantly can find 
its solution in operator form. 

At first sight the chronological operator exponent is a very specific case 
of linear operators. The following proposition asserts the opposite. Let we 
have arbitrary but linear differentiable reversible operator A{t) . From obvious 
identity 

^ = [^A-\t)]A(t), A(t)\ t=a = A{a) 

we obtain 

A(t) = f cxp{£ dr [^-A-^r^Aia) , 

so the chronological operator exponent can represent any linear differentiable 
reversible operator, only this form is exactly adjusted for solving defined DE. 
For any operator of such type there is corresponding linear DE. We can stress 
here that to solve a DE in explicit form we have to express explicitly the action 
of corresponding linear operator. It is the challenge: far from all linear operators 
can be specified by explicit action. 

In present paper for simplicity we restrict consideration by first-order DE's. 

3 Operator solutions for first-order non-linear 
ODE 

First-order differential equation of arbitrary form 

^ =/(*,«(*)) (5) 
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on account of non-linearity is not suited immediately for solving by operator 
method. Nevertheless there are many possibilities to convert the problem to 
linear one. Most of them are concerned with introduction of spaces of larger 
dimensions. Here we demonstrate some of them, which are almost generally 
applicable. 

If we introduce a new function 

S(t,w) = e uu ^ , 

from (5) one can easily derive the following equivalent to (5) first-order (with 
respect to t) linear DE for S(t,u)) 

= w /(t, A)5(t, U ) , S(a, ui) = , (6) 
which formal solution follows immediately from (3) in the form 

S(t,u)=fexp{J*dTuf(T,^ ;; )}^ a \ (7) 
The solution of equation (5) is obvious from (7): 

u{t) = y i \ u=0 . 

UU) 

We have to note here that when we look for general or any particular solu- 
tions we can utilize indefinite integrals in operator solutions. 

The operator solutions of type (3), (7) are formal enough. Meantime, solu- 
tions satisfy the following requirements to be the exact solution: 

a) . Unknown is on the left side. The right side consists of only and only 
knowns (variables, coefficients, functions, operators and etc.). 

b) . Direct substitution of the solution reduces the equation to an identity. 
The theorem for uniqueness of solution of differential equation shows that if 

we have an exact solution of a problem in several forms, all the forms can be 
transformed to each other. So the main problem that remains in this approach is 
how to transform of the operator solution to that form which can be considered 
as usable. 

The most part of this paper is devoted to hunting such ways. 



4 An idea of solving algorithm 

It would be rather well if we can find a universal algorithm, which would solve 
any DE. The most attractive feature of the operator method is unusual easiness 
in obtaining of formal (operator) solutions of DE's. Can we find on this way a 
general algorithm to obtain a standard (non-operator) solutions of DEs? 

Let us consider an idea of such algorithm on special example of solution of 
equation for unknown function u(t): 

^*) =ti » (t) + ^_ » (t)| (8) 
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where a(t) is some known function. Such type of equation is very suitable for 
testing algorithms. First of all since its particular solution is obviously known: 
u(t) = a(t). In the second place since existent CA systems do not solve at 
present most cases of such type of equations (with specific expressions for a(t)). 
The operator solution of equation (8) in form (7) is 

u(t) = f cxp{ I dtu[^ + ^ - a 2 (t)]K c Uo, (9) 

were c is an arbitrary constant (farther on in this section for simplicity let c = 0, 
i.e. we will seek here only a particular solution of equation). If we represent the 
chronological exponent as the series (4) : 

d °° 

U{t) = 5ZS*(w,t)]u;=0, 
i=l 

terms of series obviously obey the following recurrence relations 

S 1 (w,t) = l, Si(u,t)= fdtu;[^ + ^--a 2 (t)}S i _ 1 (uj,t). 



l duj 2 dt 

Now let us use the following trick. If we compute the sequence of 

U n {t) = [t— ^Si(w,f)]a, = o , 



wc find 



Ul (t) = a(t) - J dta 2 (t) , 
u 3 (t) = a(t) -2 dta(t) I dta 2 {t) + I dt[ I dta 2 {t)f 



u 5 (t)=a(t)-3 J dt[j dta 2 (t)} 2 +6 j dt j dta{t)[j dta 2 (t)} 



2 + ... 



If compare now options in sequence of u n (t) we can find that from some n 
an unaltered (survived) part in u n {t) become apparent (the a(t) here). It is 
reasonable that we can suppose that it is just the sought solution. On solution 
of DEs we always have the opportunity to verify it. If we get affirmative answer, 
so it is the end of procedure, else we proceed with larger n. 

We have prepared a computer algebra prototype of some variants of the 
algorithm using the Maple system. The algorithm manages, e.g. the following 
DE's (which we choose mostly by chance): 

DEI := =u 2 (t) + -4jt +41552420 1 90 + 4410 t S9 + - 13430 1 23 + 

+ 83 30 1 22 + 6 241 t im - 7742 1 181 - 2401 1 180 + 1, 
Ansl :=u(t) = + 79 1 91 + 49 1 90 ; 



G 



DE2 := 



du(t) 
dt 



u 2 {t) -pm cos(i) sin^" 1 ^) + p 2 sin 2m (t) 
Ans2 :=u(t) = p sin m (i); 



L>£3 := 



du{t) 
dt 




Ans3 := u{t) = pt rn e nt + q ln'(i). 



Let us make some remarks here. For specific equation we can transform 
operator solution to improve, e.g. rate of convergence. In the cases under 
consideration in this section it is very useful with help of operator identities 
(see below (10) (12)) to modify the operator solution to the following form: 



Experimentation with prototype of the algorithm leads us to conclusion that 
the algorithm is workable and can solve (at least in sense of particular solu- 
tions) a set of equations which Maple cannot solve now. It is really the solving 
algorithm as against from algorithms based on certain solvability conditions or 
mapping into separable equations. 

Central failure of this algorithm consists in the fact that it heavily relies on 
basic CAS procedures and requires noticeable resources. 

It is clear that we can try to find some particular solutions (under different 
values of c) or even general solution of DE by the same way and exploit any 
form of operator solutions. However the summation of functional series with 
several variables is for existent CAS substantially more resources consumimg 
task. 

In next sections we examine some ways for solution transformations, which 
would allow to simplify procedure of solutions calculation. 

5 Basic transformation identities 

Our aim is to reduce the equation solution from operator form (3) to more 
practical, calculable or usual (standard) expression. To do so it is necessary to 
find ways to evaluate chronological operator exponents or convert it to combi- 
nations of some operators with relatively simple properties, which would permit 
to carry out all the operators and obtain as a result the solution of equation 
" by quadratures" . 




da(t) 

dt 
da(t) 

dt 
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One of the main tool kits in the operator method is a family of identities for 
chronological operator exponents for solutions transformations. 

At present time the most- used method to "disentangle" of operator ex- 
ponents is well-known Campbell-Bakcr-Hausdorff (CBH) expansion, which for 
chronological exponents has the following version: 

f exp{- [ T dU(mB(T)fexp{[ T dtA(0=B(r)+ f dr, [B(r), i(n)] + 

J a J a J a 

...+ Tdn r ' dn-- - r 1 dT n [[...[[B(T),A(T 1 )},A(T 2 )}...},A(r n )} + 

J a J a J a 

+ dn dr 2 --- / dr n+1 f exp{- / d£_A(£_)} x 

J a J a J a J a 

x [[. . . [[S(r),4(Ti)],ii(75)] . . .], i(r n+1 )] f cxp{ d£A(Q}, 

J a 

where [B(t), A(t\)] = B{t)A{t\) — A{t\)B{t) is the standard notation of oper- 
ators commutator. 

For purpose of this paper the following related identities are of fundamental 
importance (all operators arc supposed to be linear): 

f cxp{ I drA(T)}f exp{ [ dTB{r)} = 

J a J a 

= fcxp{[ dr [A(t ) + f cxp{ j dU(0}B(r)f Q cxp{- f d£A(0}]} 

J a J a J a 

(10) 

or in more general form 

A{t)f exp{ j dTB(r)} = 

= f exp{ jf' dr [(^1)A- 1 (t) + A(r) B(r) A~\r) }A(a) (11) 

end 

f exp{ f dT[A(T)+B(T)]}=f exp{ f dr A(t)} x 

J a J a 

x f exp{ f drfo exp{- f d£ A(£)} B(t)T exp{ f d?A(£)}}. (12) 

J a J a J a 

Correctness of these identities is easily verified by differentiations (the left and 
right sides satisfy the same differential equations with the same initial condi- 
tions). 

From identity (12) follows the identity for chronological exponents with re- 
verse ordering: 
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T cxp{ / cIt[A(t)+B(t)}} 



= fa cxp{ f drfo exp{ f d£ i(£)} B(r)f exp{- f d£i(0}} x 
x f cxp{ / dri(r)}. (13) 



We will need later an variant of identity of type (11), when we do not demand 
from stationary operator D to be reversible: 



DTexp{[ dTB(r)D} = t exp{ f dr D B(t) }D. (14) 

J a J a 

The identity (10) imparts group properties to chronological operators. We 
do not evolve here this interesting matter. 

6 Generalized shift operators 

To transform formal solutions to ordinary expressions we have to express explic- 
itly the action of corresponding linear operators. Every solved equation gives 
us an example of definite action of given linear operator. So the list of "good" 
operators is nonempty. Our goal here is to study the properties of certain class 
of operators that we will name as generalized shift operators. 

The Taylor expansion for sufficiently arbitrary function $(x) can be ex- 
pressed as 

°° 1 d k d 

k=0 

Here a can be either constant or function of arguments, which are different from 
x (in the second case ^ is replaced by ^ . If we read this expression from 
the right side, we can obtain the result of action of exponential form of shift 
operator exp{a^} on function &(x) : 

cxp {a— }®(x) = $(cxp {a-^-}x) — <I>(x + a) 

or 

d d 
exp{a(t) — }$(jc) == $(exp{a(i) — }x) = $(x + a(t)). (15) 

Let us consider shift operators with new variable y — tp(x) . It is obvious, 
that 

cxp{a(t) — }<t>(y) = exp{a(t)£^y}$Ma;)) = 
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When ip{x) is fixed function, the set of shift operators with different parameters 
a(t) forms Abelian group. Shift operators with different ip(x) do not commute. 

It is quite obvious that all examined shift operators correspond to DEs with 
separable variables. 

Let us consider now one important property of the operators of the following 
type 



f exp{^ dTf{r,x)-^}, 

which we will call as generalized shift operators. 

Let g(t, x) and f(t, x) are ordinary functions (linear operators of multiplica- 
tion by value of function). Since the commutator 

lg {t ,x),f(r,x)^} = -f(r,x)(^^) 

is ordinary function, so from CBH - expansion we can conclude that following 
expression 



/"* d - f l d 

G(t, x) = Texp{J dr /(r, x)—} g(t, x) T exp{- J dr f(r,x) — } 

is ordinary function too and we can compute it simply as 

G(t,x) = Tcxp{J^ drf(r,x)^}g(t,x). (16) 



We find that 
ft 



r d 

Tcxp{ / drf(T,x) — }g 1 (t,x)g 2 (t,x) = 

f* d - f l d 

= Tcxp{ dr f(T,x) — }g 1 (t,x)T exp{- drf(T,x) — }x 

r* d 

xTexp{y dr f(T,x)—} g 2 {t,x) = 

= [fex V {J dTf(T,x)^} 9l (t,x)}[fexp{J dr f(r,x)^} g 2 (t,x)} = 

= G 1 (t,x)G 2 (t,x). (17) 

If now expand g{t,x) into the Taylor power series with respect to x, use (17) 
and then collect series backwards, we obtain analogous to (15) very important 
property (homomorphy) : 
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G(t,x)=fexp{J drf(T,x)^}g(t,x)=g(t,fexp{J dr f(r,x)^} x). 

(18) 

Now we are ready to consider another variants of operator solutions for equation 
(5) and farther modifications of solving algorithm. 

7 Solution for first-order ODE in form of gener- 
alized shift operator 

Let us start from operator solution of equation (5) in form (7) but let us rewrite 
it for convenience in following notation 

S(t, c,w)=t cxp{ jf ' drcof(r, A)} e - c , 

where c = u{t)\ t=a - With help of identity (11) we find 

S(t,c,u>) = e wc f exp{ / dTwe^7(r,^)e uc } • 1. 

CHB-expantion gives 

e-f(r, ^ = f(r, ±) + c[/(r, £),.] + ^[[f(r, ^)MM + - (19) 
Since it is easily proved by mathematical induction that 

Qui — 1 

then from expantion of f(x) into power series and backwards summation it 
follows that 

therefore we can conclude with taking into account (19) and property of the 
shift operator (15) that 



/(T ' c+ ^ ) = cxp{ l;l }/(T ' c)cxp{ -|;l } - 



ii 



Then 



By expanding the chronological exponent with help of identity (12) and prop- 
erties (16), (17) of generalized shift operators we obtain that 



S(t,c,u) = exp{ A ^ }T cxp{ - J* dr f(r, c)j-}x 
x exp{— / dr g(r, c)} exp{u; / drG(r,c)}, 

J a J a 



(20) 



where 



and 



g(T,c)=f exp{^ dU&c)^} 



d,df(T, C ) 



dc 



G(r,c)=f cxp{J\u(tc)-^- c }f(r,c). 



If we differentiate the last expression with respect to u we find that 
dS(t, c, u) 



x ex 



S(t, c, u) + e" c exp{ A cxp{- jf dr /(r, c) ^} x 

p{— / dr g(r, c)} exp{w / drG(r, c)} / drG(T,c). 

J a J a J a 



If in the second item to recover initial operator form we receive 
dS(t, c, w) 



c 



S(t,c,u) + S(t,c,u) I drG{T,c), 

J a 



du 

inasmuch as S(t,c,uj)\ w= o = 1, then 

., . dS(t, c, uo) /'*,_, . 

u(t,c) = — 1^=0 = + J dr G[r, c), 

whence it follows that 

u(t,c)=c + J rfrf exp{^ ^ /(£, c) ^} /(r, c) 

or finally 

u(t,c) = T exp{jf dT/(T,c)A}c. (21) 
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All above consideration is nonrigourous enough. So it is very important to 
verify our conclusion by direct substitution of obtained solution into equation 
(5). After differentiation of (21) and with property (18) we can be sure that 
(21) is really the solution of the equation (5). 

So we have obtained solution of first-order equation in more simple operator 
form. We can prepare now a prototype of solving algorithm based on solution 
(21) by the way similar to considered in Section 4. It would be easily found 
that in this case algorithm can handle more complicated equations and with 
more efficiency. But it remains here (in sense of present CAS abilities) the 
complication in obtaining general solutions too. In following Section we consider 
possibilities to obtain general solution out of known particular solution. 

Let us finish this Section with brief remark about interconnection between 
first-order ODE and linear first-order PDE. 

If we introduce the function 

z(t,c)=f cM-J a dTf(T,c)£-}c (22) 
it is obvious that it satisfies equation 

If we act now by operator 

/"* d 
T exp{ J dTf(T,c) — } 

on both sides of (22), we get 

z(t, u(t, c)) = c 

- the well-known equation expressing interconnection between solutions of first- 
order ODE and linear first-order PDE. 

Besides that that we obtain (operator) solutions of the equations, a nuance 
here consists in fact that conventional derivation of this interconnection is car- 
ried out by means of geometrical considerations (method of characteristics). 
We use only algebraic transformations and this approach is applicable for more 
complicated problems. We have the chance to ascertain it in the next Section. 



8 General solution in terms of one particular 
solution 

If we are able to find some particular solutions of equation (5), then arises the 
question: can we express the general solution of the equation through these 
particular solutions? One of the possible answers is the famous solution of this 
problem by Lie [9]. In terms of operator method Lie's result is: the general solu- 
tion can be expressed by relatively simple function of set of particular solutions 
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iff operator solution (7) can be expressed as a product of simple shift operators 
of type (15). On the other hand it is well-known that, e.g., if wc know only 
one particular solution for Riccati equation the general solution can be found 
through quadratures. 

Let us analyse this question from operator method point of view. 

From solution (21) and with help of identities (11), (14) we find that 



du(t, c) 
dc 



f exp{^ dT^f(T,c)}-l = 

f exp{ jf* dr [ /(r, c) + f(r, c)^] } • 1 = 

eM^drfo exp{ J\f^,c)^}^^-}. 



If denote 



then on account of (21) we can find that u(t, c) satisfies the following integro- 
diffcrential equation: 



du(t, c) 
dc 



= cxp{ / dr $(r, u(r, c)) }. 

J a 



(23) 



Solution of such type of equations is the challenge. 
Let us introduce the following functional 



/oo 
d£o;(Ou(£,c)}. 
-oo 



Since 



/OO 
-oo 



dc ' J_ 00 dc 

then to the integro-diffcrential equation (23) corresponds linear equation for 
S[c, (*>{£)] with variational derivatives 

dS[c-MZ)] r °° r£ 
dc 

with initial condition 



r°° n s 

= J ^4)exp{j d T $( T _)}S[c; W (e)] 



/OO 
d£u(£)u (£)}, uo(t) = u(t,c)\ c = c 
-OO 



Its operator solution is 
S[c;lo(0} = 



f>c f-co r-c, r /■ oc 

= f c cxp{ / da / dew(0 cxp{ / dr$(T,— -)}}exp{ / dew(0«o(0>- 
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It is obvious that solution in this form is very intricate. So we have to 
simplify it. If we repeat the chain of transformations almost literally like that 
in preceding Section, we obtain the solution in more simple form 



u(t,c) = exp{(c-c ) / d£ cxp{ / dr $(r, u (r)) }- — —-} u (t). (24) 



Now we have to examine that (24) is really the solution of equations (23) and 
(5). The verification is not difficult matter here. 

So we conclude that if given first-order ODE has general solution it can be 
expressed in terms of one particular solution of the equation. Certainly we can 
find the well-known result for Riccati equation from (24). 

The expression (24) represents the algorithm for computing the general so- 
lution from known particular solution. We have not tried to experiment with 
this algorithm yet. 

9 General solution in terms of arbitrary func- 
tion and solvability conditions 

Let now f(t, c) — g(t, c) + h(t, c), then we can find the following form of general 
solution of equation (5) 



u(t, c) = T exp{ / drf(T,c) —} c = f cxp{ / dr [g(r, c) + h(r,c)] —} c = 






5 




= f exp{ f dr [f cxp{ f d£ h(£, c) |-} g(r, c)} x 



J a J a 




or 



u 



(t, c) - T cxp{ f dr [To exp{ f d£ h(£, c) |-} [/(r, c) - h(r, c)}] 



J a J a 



X 




(25) 



where we can consider h(t, c) as an arbitrary differentiable function. 
If it is known the general solution 



z(t, c) = T exp{ / drh{T,c) — }c, 
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then with help of property (18) we obtain that 



u(t, c) = T exp{ / dr 




dz{r,c) 



g- c }z(t,c). 



(26) 



As far as h(t, c) is an arbitrary function, therefore preceding expression is valid 
for any differentiable function z(t,c). The particular solution of equation (5) 
with initial condition u(t, c)\ t=a — c is obtained when z(t, c)\ t=a — c. 

In spite of obviousness and simplicity the form (26) has important meaning. 
And first of all it gives some recipes for finding ordinary (standard) solution of 
the initial equation (5). 

In contrast to all above mentioned operator forms of solution the last de- 
scribed is not algorithmical as long as z{t, c) remains to be an arbitrary function. 
We do not analyse completely here all consequences of interpretations of solu- 
tion in form (26). Let restrict ourselves by consideration of some solvability 
conditions only. 

Having noted that if 



corresponds to right-hand part of solvable first-order ODE, then the general 
solution of given equation can be obtained too. 

If we require $(i, c) = 0, we need to solve the same given equation. In 
another cases we will need to solve different equation and have chance really 
to solve it. Main difficulties on such a way arise in the fact that we can not 
circumscribe all solved DE by single expression, so we have the dilemma: if we 
choose a certain z(t, c) we are not sure that the equation with corresponding 
$(i, c) is solvable and vice versa - we are not sure that the equation (27) in 
respect to z(t, c) with some <!>(£, c) selected in advance is solvable. 

The situation becomes simpler if we content ourselves with some solvability 
conditions. It is clear that we are interested to pose the broader solvability 
conditions. 

To formulate an algorithmical solvability conditions we have to fix some 
finite set of functions z(t, c) - the set of trial functions. Then by using some 
criterion we successively examine these function trying to find at least one of 
them (generally speaking, there are a family of such functions) that satisfy the 
criterion. If such a function is found, then solution of given equation is obtained 
almost automatically by substitution of this function to (26). 

The following important classes of equations that permit obtaining of general 
solutions as a whole can illustrate this approach. 

Linear equations. Here $(i, c) = a(t)c + [3(t). It is obvious that solvability 
condition in this case can be formulated as follows. If z(t, c) is found to be 
satisfy the criterion 



*(*,c) 



f{t,z{t,c))- 



dz(t,c) 
dt 



(27) 



dz(t,c) 
dc 




dz(t,c) 
dt 



} = o, 



(28) 



Or 
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then the given equation can be solved through the function z(t, c) and a linear 
equation. Since in this case u(t,c) — A(t)z(t,c) + B(t) then here we use linear 
transformations of variable u(t, c). Such type of solving method is due to Liou- 
ville. Relatively recent implementations of this approach to CAS make a good 
showing in solving many of Kamke examples [10]. 

Equations with separable variables. Here $>(t, c) = a{t)j3(c). The solvability 
condition in this case can be formulated as follows. If z(t, c) is found to be 
satisfy the criterion 

82 M ^L-^ H-o. <») 



dtdc L L gg&g) 

then the given equation can be solved through the function z(t, c) and an equa- 
tion with separable variables. The suitable functions a(t) and /3(c) under known 
z(t, c) are easily obtained from (27). 

Reduction of an equation to the one with separable variables of course is 
the well-known method, but as a matter of fact for solving non-linear equations 
it is used mainly heuristically simply by guessing appropriate substitution. In 
resembling a way L. Chan and E.S. Cheb-Terrab have implemented recently so- 
colled hypergeometric solutions method [11] where the sets of functions z(t,c) 
and substitutions are fixed and only second-order linear (or Riccati) equations 
are considered. In contrast to their method we here have to fix only the set of 
trial functions z(t, c) so our approach can handle broader family of equations. 

Bernoulli equations. Here $(t, c) = a(t)c v + [i(t)c (i] ^ 1). The solvability 
condition in this case can be formulated as follows. If z(t, c) is found to be 
satisfy the criterion 

a r /( M(t , c ))-g%£) 1 _ nt ,z(t,c))-S^l 

i-{ =^= ± } = 0, (30) 

dc c"! 

then the given equation can be solved through the function z(t, c) and an 
Bernoulli equation. The suitable functions a(t) and (3{t) under known z(t, c) 
are obtained from (27). 

Despite the fact that transformations here are simple enough we do not know 
any real solving method based on this scheme. 

We can easily find examples of equations that are solved in such a way and 
are not solved say by Maple. 

It is obvious that list (28)- (30) is far from to be complete, analogous solv- 
ability conditions can be obtained from another solved equations. And what 
is more the solution (26) can be iterated what lead to schemes with successive 
variables changing (substitutions). In other words we can use combinations of 
conditions too. The general strategy for such type of algorithms needs further 
development. 



17 



10 Conclusions 



We presented some ways in solving DE by operator method. Besides first- 
order DE we can obtain operator solutions for linear and non-linear equations 
of arbitrary order. The properties of its solutions expect further studies. 

In conclusion we believe that we succeeded in demonstration the fact that 
operator forms of DE solutions can be handled analytically none the worse 
than ordinary functions. In some cases its transformation properties are more 
comfortable than, e.g., for some special functions. 
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